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We find that the probability distribution for the 
for different systems including random walk and 
infinite set of singularities at I = 1/k with k = 2, 
additionally, p(l) has an essential singularity at 
dimensional situation. 
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largest intervals p(l) exhibits universal properties 
random cutting models. In particular, p(l) has an 
3, . . . which become weaker and weaker as k — > oo; 

1 = 0. Phese properties are found also in many 



Systems with complex landscapes are of great current 
interest in a variety of fields ranging from optimization to 
^f^evolution [1]. The phase space in such systems is usually 



'broken into many valleys. The multivalley structure de- 
termines many peculiar properties as it was first realized 
for spin glasses [2] . The large valleys fluctuate from sam- 
^ple to sample and lead to anomalous long-time behavior, 
' 'ergodicity breaking, etc. Large valley distribution and 
^Oother extremal characteristics play therefore a significant 

role in various problems. 
ON To study extremal properties of random systems, we 
^Pconsider two very simple geometrical models. The first 
j~ model arose as a byproduct of our previous work [3] which 
f~^studies the kinetics of annihilation of charged particles in 
two-component plasma with a truly ID Coulomb in- 
^^teraction and neutrality condition. This system is decou- 
led into a succession of uncharged chains which evolve 
ndependently [3]. The long time behavior is therefore 



^determined by long chains, in particular the time of the 
F^complete annihilation is directly related to the longest un- 
r-|charged chain. Thus we run into a problem of the longest 
Ouncharged chain distribution which exhibits a surpris- 
ingly rich behavior. Although some of the features are 
model-dependent, the basic properties of the longest in- 
terval distribution function such as the existence of an 
infinite countable set of singularities and their locations 
seem universal. At least, we will demonstrate that simi- 
lar features arise in another very different model, namely 
the random cutting model. We therefore anticipate that 
our qualitative conclusions should provide an insight on 
the probability distribution of the largest valley in disor- 
dered systems and systems with a complex landscape. 

In the first problem, we consider sequences of 2N 
charges Qj, Qj = ±1, with ^2Qj = to satisfy the 
neutrality condition. It is clear that such a sequence can 
be visualized as a random walk (RW) starting at the ori- 
gin at time t = 0, jumping a step up at t = j if Qj = 1 
and a step down if = —1, and returning back to the 
origin at t = 2N (Fig. 1). A RW can return to the origin 
few more times, and we ask for the probability distribu- 



tion of the longest time interval, 2L, between successive 
departures and arrivals (see Fig. 1). 




FIG. 1. A typical random walk with 2N = 100 and 2L = 46 



A typical number of returns scales as y N , and so does 
a typical time interval. However, the longest time interval 
scales as N rather than \/N. To demonstrate this and 
other peculiar properties of the longest time interval let 
us calculate the total number, T^(L), of RWs with the 
longest time 2L. If 2L > N , a time interval of length 2L 
is always the longest one and therefore T^(L) is easily 
computed: 
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Here Rk = (2K)\/ K\K\ gives the total number of "head" 
RWs of length 2K from the origin (t,x) = (0,0) to 
the beginning of the longest interval (2K, 0), Fl = 
2{2L - 2)1/ (L - l)\L\ gives the total number of RWs 
starting at {2K, 0) and then first returning to x = 
at t = 2K + 2L, and i?jv-if-i gives the total number 
of "tail" RWs of length 2N - 2K - 2L [4]. Computing 
the sum in Eq. (1) yields T N {L) = F L ■ 2 2N ~ 2L and thus 
the longest interval probability, P^(L) = Tn(L)/Rn, is 
equal to F L ■ 2 2N ~ 2L /R N for 2L > N. In the continuum 
limit — 5- cxd it is convenient to use the probability den- 
sity p(l) = NPn(L) where / = L/N is the scaled length. 
Our previous result can then be rewritten as 
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When 2L < N , we should modify Eq. (1) to guarantee 
that the head and tail regions do not contain intervals 
of length > 2L. In the continuum limit, the probability 
distribution p(l) obeys the integral relation: 
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Here x = K/N and 1 — / — x are the scaled lengths of the 
head and tail regions, /(/, x) and /(/, 1 — / — x) are the 
probabilities that the longest intervals in the head and 
tail regions are smaller than /. For example, 
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when x > / while for x < I we have /(/, x) = 1. Besides /- 
subf actors, the first and the last factors inside the integral 
in Eq. (3) are just the continuum representation of Rk 
and i?jv-if-i, respectively. 

Notice that /(/, x) depends on the the tail of the prob- 
ability distribution p(rj) for rj > l/x > 1/(1 — /). In 
particular, when I >\, /-subfactors are trivial, /(/, x) = 
/(/, 1—l—x) = 1 and Eq. (3) reduces to Eq. (2). Similarly, 
in the interval i < / < i we readily compute /-subfactors 
from the knowledge of the probability distribution in the 
preceding interval, Eq. (2). After elementary computa- 
tions we find 
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for i < / < 
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These results can in principle be generalized for any 



< 



k + l 

cumbersome 



< j but even for k 



3 the outcome is rather 
Let us just stress the most important 
feature: Since the expressions for the probability dis- 
tribution function p(l) are different in different regions 
^q-j- < / < i, singularities are expected at / = l/k. These 
singularities, however, become weaker as k increases. Ex- 



panding expressions (2) and (5) near . 



finds that 



they differ by a term p s ing(0 ~ |/ — 5 1 3 / 2 - Similarly, ex- 
panding (5) and p(l) in the next region 5 < / < g near 
/ = i one gets p s ing(0 ~ |/ — ^| 3 . The nature of singu- 
larities then becomes clear, 



Psing (0 
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Even the first singularity at / = -| is hardly visible (see 
Fig. 2). However, numerical differentiation of the data 
allows us to detect it. 



FIG. 2. Continuous probability distribution p(l) for the 
RW problem. Numerical simulations of 10 8 different config- 
urations of length 2N = 1000. Values of p(l) are given by 
Eq.(2) for / > 1/2 and Eq.(5) for 1/3 </ < 1/2. 

Our numerical results also indicate that the probability 
density has an essential singularity p(l) ~ exp(— const//) 
when / — y 0. To explain this behavior heuristically note 
that if the longest time interval between the successive re- 
turns of a RW to the origin is 2L, this RW typically stays 
within a region (—\/T, \/T). The probability of that can 
be estimated by solving the diffusion equation in a cage 
with adsorbing boundary conditions. In a cage of size B 
this survival probability behaves as exp(—t/B 2 ). In our 
case t = N, B 2 ~ L and thus p(l) ~ exp(— const//). 

Notice that a closely related problem has been re- 
cently investigated [5]. In that model, the maximal- 
length chains were different from ours (just uncharged, 
not necessarily starting at the origin, and making arbi- 
trary number of returns to the starting point). Numerical 
simulations [5] demonstrate a singularity at / = i, even 
more sharp than ours. The next singularity at / = i 
has not been notified. However, all the singularities do 
exist and are located at the very same places, / = i, as 
in our model. The singularities just become weaker and 

weaker (p s ing(0 ~ \l — j\ seems possible) and hence 
are hardly visible for k > 3. 

We now turn to the random cutting model which is 
defined as follows. Starting with the unit interval we 
cut it (the cutting probability is assumed to be uniform). 
Two intervals are formed, left and right, and in the next 
stage we cut the right one only. The cutting process is 
then continued with the right interval from the preced- 
ing stage. We again wish to determine the probability 
density, p(l), for the longest interval which is formed in 
the cutting process. 

When / > i, the probability density satisfies an inte- 
gral equation 



p(l) = 1 + 
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Indeed, we can get the desired interval of length / at the 
first stage of the process (since I > | we never get a 
piece of length / or larger at any further stage). This 
gives the first constant contribution to the right-hand 
side of Eq. (7). A complimentary possibility means that 
at the first stage we cut at some point x with x < 1 — 
/. The probability density to get the desired interval of 
length / from remaining interval of length 1 — x is given 
by jzr^p (t^?) > where the prefactor provides a proper 
normalization. 

In general, the governing integral equation becomes 

f 1 -' ( x \ dx [' ( I \ dx 

(8) 

The last integral term is built as previously up to a small 
correction - the upper limit is now / to guarantee that 
the first cutting gives the interval x less than /. The two 
first terms again describe the situation when we get the 
/-length interval at the first stage and guarantee that in 
the following cutting of the remaining interval of length 
1 — /, we never get an interval of length / or larger. For 
/ > \, Eq. (8) is identical to Eq. (7) if we agree to put 
p(x) = for x > 1. 

Changing variables inside the integrals, y = y^y in the 
former and z = in the latter, simplifies Eq. (8) to 

f 1 dz 
p(l) = 1- / P(y)dy+ / p(z) — . (9) 

Differentiating Eq. (9) we finally reduce it to the 
difference-differential equation: 

When / > i we have dp/ dl = —p/l with the boundary 
condition p(l) = 1. This is readily solved to find 

p(l) = j, for l>\. (11) 

If i < / < i, the variable 1/(1 — I) belongs to the region 

[■|,1]. Thus p (^jzr^j = ^~r~> anc ^ (^) becomes a 
simple differential equation which is solved to find 

MO^ + yln^), for \<l<\. (12) 

It is possible to derive recursively exact expressions for 
the probability density in consecutive regions < / < 
The main conclusion remains the same: There are 
singularities at / = ^ which become weaker and weaker 
as k increases. More precisely, 
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FIG. 3. Probability distribution p(l) for the cutting process 
in d = 1, 2, 3. Numerical simulations of 10 9 processes. 

Consider now the small-/ behavior. Our simulation 
results indicate an essential singularity at / = 0. Let us 
try the same type of singularity, p(l) ~ exp(— a/l), as 
in the first problem. Substituting it into Eq. (10) gives 
different power law pref actors, l~ 2 in the left-hand side 
and / _1 in the right-hand side. This defect can be readily 
removed by replacing the constant a by a slowly varying 
function of /. An asymptotic analysis yields 

m ( ln(l//) + lnln(l//) + .. A 
p(/)~expf K -L-L j^J^ . (14) 

Our simulational results well agree with the asymptotic 
prediction of Eq. (14). 

The singularities seem to be a generic feature of the 
random cutting model and the likes. It may be diffi- 
cult to prove generally, but we believe that the mere 
existence of singularities and their location are robust 
to small changes in the model rules. In particular, if 
one considers not necessarily uniform cutting probability, 
p(x) ^ const, one finds that the present approach still ap- 
plies and only quantitative changes appear, e.g., integral 
equations will contain the cutting probability p(x). To 
verify the robustness of generic properties we consider a 
(/-dimensional random cutting model. For example, in 
two dimensions the model is defined as follows. We start 
with a unit square, pick an arbitrary internal point and 
cut the square along the two lines going through this 
point in x- and ^-directions. We then continue the cut- 
ting process with only one of four rectangles, say with 
the north-east one. We wish to determine the probabil- 
ity density, p(A), for the rectangle of the largest area A 
which arises in the cutting process. 
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Note first that for the uniform cutting probability, 
p(x, y) = 1, the area distribution function is not uniform 
even in the one-cut process. Instead, it is given by 

p {A) = f f dxdyS(xy-A) = -ln(A). (15) 
Jo Jo 

The probability density satisfies an integral equation 

p(A) = 3 Po {A) + [ [ dxdy 6(xy — A)— p (— \ , 
Jo Jo X V \ x vJ 

(16) 

where 6 is the step function and the factor 3 in the 
first term accounts for three rectangles which are kept 
unchanged in the following cutting process. Eq. (16) 
is derived in the same way as Eq. (7) in the corre- 
sponding ID version and valid for A > again like 
in the ID situation. Changing variables from (x,y) to 
(u = xy,v = x/y), performing the ^-integration, twice 
differentiating the remaining integrals, and solving the 
resulting differential equation, we find an exact solution, 

P( A ) = l(j- A )> fOT 

This result perfectly agrees with simulations (Fig. 3). 

Simulations in two dimensions also clearly indicate sin- 
gularities at A = i and j. However, we still expect the 
singularities at all inverse integers, A = ^. The qualita- 
tive origin of singularities is the same as previously. If in 
the cutting process a sufficiently large rectangle, namely 
a rectangle of area A > \ is formed, it is always the 
largest. For the largest rectangle of a smaller area A < \ 
we should additionally guarantee that other rectangles 
are smaller than A and hence the probability distribu- 
tion function p(A) is described by different expressions 
for A > 77 and A < \. When we pass through other crit- 
ical points A = ^ new possibilities arise and hence the 
probability distribution function changes again. Geome- 
try comes into play and changes the strength of singular- 
ities, but neither their existence nor location. Indeed, in 
2D one can form two rectangles of the area A = \ or four 
rectangles of the area A = \ just in one cut. However, 
one cannot form three rectangles of the area A = ^ after 
one cut - one needs an additional cut. Note that a singu- 
larity which arises at the later stage of the cutting process 
is weaker - mathematically it appears by "integrating" 
the previous singularity. In one dimension the strengths 
of singularities were consecutive, see Eqs. (6) and (13), 
since they appear on the consecutive stages of the cutting 
process. In 2D both singularities at A = \ and i arise 
after the first stage and thence identical (simulations in- 
dicate that p smg (A) ~ \A - \\ and p smg (A) ~ \A - \\). 
The singularity at A = ^ appears later and should be 
weaker. Unfortunately, we could not find its strength 
neither theoretically nor numerically. 



We believe that the largest volume distribution shows 
the same properties in all dimensions. We have veri- 
fied this in three dimensions. The singularities become 
weaker as the dimension grows but we have certainly de- 
tected first few of them. 

In summary, we have demonstrated surprisingly pe- 
culiar properties of the longest-intervals probability dis- 
tribution. We have found an infinite set of singularities 
at the inverse integer values of the the longest interval 
/ = i. Singularities become weaker and weaker as k in- 
creases. Additionally, we have found an essential singu- 
larity at / = 0. We have confirmed the existence of these 
singularities for several models. For the random cutting 
model, these results are valid in arbitrary dimension. We 
expect these singularities to appear in other apparently 
unrelated problems. 
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